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ABSTRACT 

We present the ellipticity distribution and its evolution for early-type galaxies in clus- 
ters from z ~ 0.8 to the current epoch, based on the Wide-field Nearby Galaxy-cluster 
Survey (WINGS) (0.04 ^ z ^ 0.07), and the ESQ Distant Cluster Survey (EDisCS) 
(0.4 ^ z ^ 0.8). We first investigate a mass limited sample and we find that, above 
a fixed mass limit (M* ^ lO^^'^M©), the ellipticity (e) distribution of early-types 
noticeably evolves with redshift. In the local Universe there are proportionally more 
galaxies with higher ellipticity, hence flatter, than in distant clusters. This evolution 
is due partly to the change of the mass distribution and mainly to the change of the 
morphological mix with z (among the early types, the fraction of ellipticals goes from 
~ 70% at high-z to ^ 40% at low-z). Analyzing separately the ellipticity distribution 
of the different morphological types, we find no evolution both for ellipticals and for 
SOs. However, for ellipticals a change with redshift in the median value of the distribu- 
tions is detected. This is due to a larger population of very round (e < 0.05) elliptical 
galaxies at low-z. In order to compare our finding to previous studies, we also assem- 
ble a magnitude- "delimited" sample that consists of early-type galaxies on the red 
sequence with —19.3 > Mb + I.2O82; > —21. Analyzing this sample, we do not recover 
exactly the same results of the mass-limited sample. his indicates that the selection 
criteria are crucial to characterize the galaxy properties: the choice of the magnitude- 
"delimited" sample implies the loss of many less massive galaxies and so it biases the 
final conclusions. Moreover, although we are adopting the same selection criteria, our 
results in the magn itude- "delimited" sample are also not in agreement with those of 
iHolden et al.l ()2009f) . This is due to the fact that our and their l ow-z samples have a 
different magnitude distribution because the iHolden et al.l ()2009l ) sample suffers from 
incompleteness at faint magnitudes. 
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1 INTRODUCTION 

Ellipticals and lenticulars (SOs) belong to the class of early- 
type galaxies. This means that they have several properties 
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in common: they dominate the total galaxy population at 
high masses, they preferentially inhabit dense regions of the 
universe, such as rich clusters (jPressler et aljriQQTl ). they 
tend to be passive, they have red colors and their spectra 
sho w strong values of the characteristic Z?4000 featu re (see 
e.g. iKauffmann et al.l [200^ : iBrinchmann et al.|[20o3 ): they 
lack spiral arms and in most cases exhibit neither major 
dust features nor a large interstellar gas content. For these 
reasons, often they are considered together. 

On the other hand, elliptical and SO galaxies differ in 
several important ways: SOs are bulge-dominated systems 
with an identifiable disk (e.g. lScorza fc van den Boschll 19981 : 
iLaurikainen et al] l2007f). that is mainly rotationall y sup- 
ported fe.g. lErwin et aLli2003l : [Cappellari et alll2005l'). their 
intrinsic shape is similar t o that of spirals I Rood fc BaumI 
1 19671 : ISandage et al.' 'l970') and their formation is still not 
well understood. [Hubble ( 1936) first proposed their exis- 
tence as transitional class between ellipticals and spirals. 
Understanding how they form and evolve is essential if we 
wish to have a complete picture of how galaxy morphology 
is related to galaxy formation and the environment. 

Then again, ellipticals show ellipsoidal shapes, not 
rarely with significant kinematic twists, and kinematically 
decoupled components in their centr es. Most of them are 
not c haracterized by strong rotation l|Bertola fc Capacciolil 
Il975f l. and their luminosity profiles follow a Sersic's law. In 
the Local Universe disky ellipticals are probably the high 
bulge mass end of SO. 

Morphologically, iDressler et al] j 19971 ) showed that, at 
least for bright galaxies, the raising fraction of early-type 
galaxies since z ~ 0.5 corresponds mainly to an increase 
of lenticular SO galaxies, with a roughly constant ellipti- 
cal fraction. SOs are quite rare in clusters at high redshift 
(z > 0.3 — 0.4); as a consequence, they have to acquire their 
shapes with different time-scales and later than ellipticals. 
The evolving fraction of SOs in clusters might result from the 
evolving population of newly accreted spiral galaxies from 
infa lling groups and the field. 

iFasano et ahl |200Cf) showed that the cluster SO to ellip- 
tical ratio is, on average, a factor of ~ 5 higher at z ~ than 
at z ~ 0.5. At higher redshift, there is no evidence for any 
further evolution of the SO fraction in cluster s to 2 '--^ 1: most 
of th e evolution occurs since z ~ 0.4 (see e.g . I Post man et al.l 
12005^: 'Pesai et al."2007': Wilman ct al. 2009 



iDressler et al.. (.1997. ') and Postman et all (|2005l ) also in- 
vestigated the ellipticity distributions of the SO and ellipti- 
cal galaxies in their magnitude limited samples. They found 
that the ellipticity distribution of SO and elliptical galaxies 
show no evolution over the broad redshift ranges in their 
samples. Moreover, they differ from each other, providing 
evidence for the existence of two distinct classes of galaxies. 

In contrast, in their magnitude- "delimited" sample 
(with both an up per and a lower magnitude limit), 
iHolden et al.l (|2009l ) found no evolution in neither the me- 
dian ellipticity nor the shape of the ellipticity distribution 
with redshift for early-type (ellipticals -I- SOs) red-sequence 
galaxies. This lead them to conclude that there has been 
little or no evolution in the overall distribution of bulge- 
to-disk ratio of early-type galaxies from z ~ 1 to 2: ~ 0. 
Assuming that the intrinsic ellipticity distribution of both 
elliptical and SO galaxies separately remains constant, they 
finally concluded that the relative fractions of ellipticals and 



SOs do not evolve from 2 ~ 1 to z; = for a red-sequence 
selected sample of galaxies. 

All the cited works analyzed samples limited in some 
ways by magnitude cuts. For the first time, in this paper we 
analyze the evolution of the ellipticity distribution of early- 
type galaxies also in a mass-limited sample. For the sample 
in the Local Universe, we analyze the data of the Wlde- 
field Nearby Galaxy-cluster Survey (WINGS) l|Fasano et al.l 
I2OO6V while for that in the distant Universe we use the ESO 
Distant Cluster Survey (EDisCS) (jWhite et al.ll2005l ). These 
large cluster samples and their high quality images (see §2) 
allow us to characterize properly the cluster environment 
at the two redshifts and to subdivide galaxies into the dif- 
ferent morphological types and obtain robust estimates of 
ellipticity. 

This paper is organized as follows: in jj2 we present the 
cluster and g alaxy samples (W INGS (|Fasano et al . 20QJ) 
and EDisCS j White et al.l 120051 ) describing the surveys, 
the data reduction, the determination of morphologies, el- 
lipticites and masses. We also depict the selection criteria 
we follow to assemble the mass-limited and the magnitude- 
"delimited" samples. In §3 we show the results of our anal- 
ysis of the evolution of the ellipticity distribution with red- 
shift in our mass-limited samples, while in §4 we show the 
same for the magnitude- "delimited" samples. In §5 we try to 
reconcile the results of the different samples, while in §6 we 
compare our results with those found in literature ( in partic- 
ular with the results drawn bv lHolden et"ai]|2009l ) . Finally, 
in §7 we discuss and summarize our findings. 

Throughout this paper, we assume Ho = 
70kms"^Mpc"\ Q.m = 0.30, = 0.70. The adopted 
initial mass function is a lKroupal (|200ll ) in the mass range 
0.1-100 Mq. 



2 CLUSTER AND GALAXY SAMPLES 

To perform the study of the ellipticity (e = 1 — b/a, b = semi- 
minor axis, a = semi-major axis) distribution and its evolu- 
tion from z ^ 0.8 to z ~ for early-type galaxies and for el- 
lipticals and SOs separately, we assemble two different galaxy 
cluster samples in two redshift intervals: we draw the sam- 
ples at low-z from t he Wide-field N earby Galaxy-cluster Sur- 
vey (WINGS) (Fasano et al.ll200di ') and those a t high-z from 
the E SQ Distant Cluster Survey (EDisCS) jWhite et all 
I2OO5I ). 

First of all, we use a mass-limited sample, that ensures 
completeness, i.e. includes all galaxies more massive than 
the limit regardless of their color or morphological type. We 
think that this is the best choice to characterize properly 
galaxy properties. 

Then, since iHolden et al] (|2009l ) have already analyzed 
the ellipticity distribution using a sample delimited in mag- 
nitude both at faint and bright magnitudes, in order to com- 
pare our results with theirs, we also assemble a magnitude- 
"delimited" sample, following their selection criteria. 
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2.1 Low-z sample: WINGS 

The main goal of WING^I jPasano et all 120061 '). a multi- 
wavelength survey of clusters at 0.04 < z < 0.07, is to char- 
acterize the photometric and spectroscopic properties of 
galaxies in nearby clusters and to describe the changes of 
these properties depending on galaxy mass and environ- 
ment. The project was ba sed on deep optica l (B, V) wide 
field images of 77 fields jVarela et al.l |2009| ) centered on 
nearby clusters of galax;ies selected from three X-ray flux 
limited samples compiled f rom ROSAT All- Sky Survey data 
l|Ebeling et al.lll998l . lioool) and the X-ra y Brightest Abell- 
type Cluster sample ( Ebeling et al.|[l996l '). 

WINGS clusters cover a wide range of velocity disper- 
sion (Tcius (typically 500 — 1100 km s~^) and a wide range of 
X-ray luminosity Lx (typically 0.2 — 5 x lO'^^erg s~^). 

The survey has been complemented by a near-infrared 
(J, K) survey of a s ubsample of 28 clusters obtained with 
WFCAMQUKIRT (IValentinuzzi et al.l |2009| ), by a spec- 
troscopic survey of a subsample of 48 clusters, obtained 
with the spectro graphs WYFFOS@WHT and 2dF@AAT 
l|Cava et al.ll2009l ). and by U broad-band and Ha narrow- 
band imaging of a subset of WINGS clusters, obtained with 
wide-field cameras at different telescopes (INT, LBT, Bok) 
ijOmizzolo et al.ll2010l '). 

The spectroscopic target selection was based on the 
WINGS B, V photometry. The aim of the target selection 
strategy was to maximize the chances of observing galaxies 
at the cluster redshift without biasing the cluster sample. 
Galaxies with a total 1/ ^ 20 magnitude, a V magnitude 
within the fiber aperture of V < 21.5 and with a color within 
a 5 kpc aperture of [B — V)5kpc ^ 1.4 were selected, to re- 
ject background galaxies. The exact cut in color was varied 
slightly from cluster to cluster in order to account for the 
redshift variation and to optimize the observational setup. 
These very loose selection limits were applied so as to avoid 
any bias in the colors of selected galaxies. 

Our optical imaging covers a 34' x34' field. This imaging 
corresponds to about 0.6 ii2oo or more, for most clusters, 
although in a few cases only ~ 0.5ii2oo is covered. -R200 
is defined as the radius delimiting a sphere with interior 
mean density 200 times the critical density of the Universe at 
that redshift, and is commonly used as an approximation for 
the cluster virial radius. The -R200 values for ou r structures 
are co mputed from the velocity dispersions by ICava et al.l 
(|2009l) . 



2.1.1 Morphologies 

Morphological types are derived from V-band images us- 
ing MORPHOT, an automatic tool for galaxy morphology, 
purposely devised in the framework of the WINGS project. 
MORPHOT was designed with the aim to reproduce as 
closely as possible visual morphological classifications. 

MORPHOT extends the classical CAS (Concen- 
tratio n/ Asymmetry /clumpinesS) parameter set jConselicel 
I2OO3I ). by using 20 image-based morphological diagnostics. 
Fourteen of them have never been used, while the remaining 
six [the CAS parameters, the Sersic index, the Gini and M20 
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Figure 1. Normalized distributions of the MORPHOT Az- 
imuthal Coefficient for the visually classified ellipticals (366 ob- 
jects, red histogram) and SO galaxies (267 objects, blue his- 
togram) of the MORPHOT calibration sample. The Azimuthal 
Coefficient measures the correlation between azimuth and pixel 
flux relative to the average flux val ue of the elliptical iso phote 
passing through the pixel itself (from lVulcani et al.ll2010bl ). 
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coefRcients (|Lotz et al.ll2004l 'l] are already present in the lit- 
erature, although in slightly different forms. An exhaustive 
description of MORPHOT will be given in a forthcoming 
paper (Fasano et al. 2010b), where also the morphological 
catalogs of the WINGS clusters will be prese nted and dis- 
cussed . Pr ovisionally, we refer the reader to iFasano et al.l 
|2003) and iFasano et all l|2010l . Appendix A therein) for an 
outlining of the logical sequence and the basic procedures 
of MORPHOT. Here we just mention that, among the 14 
newly devised diagnostics, the most effective one in order 
to disentangle ellipticals from SO galaxies turned out to be 
an Azimuthal coefflcient, measuring the correlation between 
azimuth and pixel flux relative to the average flux value of 
the elliptical i sophote passing through the pixel itself. From 
IVulcani et al.l ( j2010b ) , we report here in Figure[T]a plot illus- 
trating the capability of the distributions of the Azimuthal 
coefficient in disentangling elliptical from SO galaxies, a cru- 
cial point in the present analysis. 

More importantly for our purposes, the quantitative dis- 
crepancy between automatic (MORPHOT) and visual clas- 
sifications turns out to be similar to the typical discrepancy 
among visual classifications given by experienced, indepen- 
dent human classifiers (r.m.s.~1.3-2.3 T types). The last 
one has been estimated from a sample of 233 SDSS galaxies 
included in the Third Reference Catalog of Bright Galaxies 
(de Vaucouleurs et al. 1991, RC3), whose visual classification 
was carried out independently by GF and AD and also com- 
pared with that given in the RC3. The comparison between 
visual (GF) and automatic (MORPHOT) classification is il- 
lustrated in Figure[5]for the MOR PHOT calibration sam ple 
(931 galaxies). In this figure ffrom Fviilcani et al.ll2010"bl ') the 
automatic classification is also shown to be bias-free in the 
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Figure 2. Comparison between visual and MORPHOT broad 
morphological types for the 931 galaxies of the MORPHOT cal- 
ibration sample. In each one of the 2D bins of the plot the per- 
centages of the visual broad types (Es, SOs, early spirals [SpE], 
late spirals [SpL] and irregulars) falling in different bins of the 
MORPHOT (broad) classification are reported on the top. Simi- 
larly, on the right side of each bin the percentages of MORPHOT 
types falling in different bins of the visual types are reported. 



overall range of morphological types, perhaps apart from 
the last bin, i.e. that relative to the very late and irregular 
galaxies. 

For now, we can apply MORPHOT just to the WINGS 
imaging, because the tool is calibrated on the WINGS imag- 
ing characteristics, and we defer to a later time a more 
generally usable version of the tool. In the following, for 
EDisCS we will use visual morphological classifications. To 
verify directly that the two methods adopted at different 
redshifts (see Section \2.2\ are consistent, we can apply the 
same "method" (visual classification and persons) that was 
used at high-z on the low-z images. 

To this aim, 3 of the classifiers that in 2007 visually 
classified all the EDisCS galaxies (BMP, A AS, VD) now 
performed a visual classification of WINGS galaxies. This 
was done on the subset of WINGS galaxies that was used to 
calibrate MORPHOT on the visual WINGS morphologies, 
including only galaxies that enter the sample we analyze in 
this paper (173 galaxies). 

The results (see Figure O from IVulcani et all l2010bl ) 
show agreement between the three broad morphological 
classes assigned by the EDisCS classifiers with the WINGS 
visual classification in ~83% of the cases, and with MOR- 
PHOT in ~75% of the cases. Again, these discrepancies 
turn out to be similar to the typical discrepancy among 
visual classifications given by experienced, independent hu- 
man classifiers, so we conclude that the different methods 
adopted provide a comparable classification. 



Figure 4. Comparison between the ellitpicity estimation from 
GASPHOT and GIM2D for early-type galaxies of the cluster 
A119. 

2.1.2 Ellipticity measurements 

Galaxy ellipticities have been computed using the tool GAS- 
PHOT (jPignateUi et al.i r2006). This is heavily based on SEx- 
tractor ("Source Extrac tor") galaxy photometry package 
(IBertin fc Arnouts|[l99^ and provides, among other quan- 
tities, ellipticity profiles of galaxies extracted from CCD 
frames. It fits simultaneously the major and minor axis light 
growth curves of galaxies with a 2D fiattened Sersic-law, 
convolved by the appropriate, space-varying point-spread- 
function (PSF), which was previously evaluated by the tool 
itself using the stars present in the frame. This approach 
exploits the robustness of the ID fitting technique, saving 
at the same time the capability, typical of 2D approaches, of 
dealing with PSF convolution of flattened galaxies. The tool 
was previously tested for non-Sersic profiles and blended ob- 
jects and its results com pared with other tools, such as GAL- 
FIT (iPeng et all |2002|) and GIM2D (iMarleau fc SimardI 
Il998l l. as shown in lPignateUi et al.1 (j200ff [ 

Since in our analysis we are comparing the ellipticities 
of WINGS galaxies to those of EDisCS galaxies, that have 
been determined using the tool GIM2D (see Section [22]), we 
also performed a comparison between the values estimated 
by the two different tools, for early-types galaxies in the 
WINGS cluster A119. As shown in Figure HI the estimates 
are in good agreement (r.m.s.~0.07). 

2.1.3 Galaxy stellar masses 

Stellar masses have been determined using the relation 
bet ween M/Lb and res t-frame {B — V) color, follow- 
ing [Beir^^^eJonJ l|200lh . who used a spectrophotometric 
model finding a strong correlation between stellar mass-to- 
light M/L ratio and optical colors of the integrated stellar 
populations for a wide range of star formation histories. This 
method was chosen to be consistent with that adopted for 
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Figure 3. Comparison between tiie visual classification performed by the EDisCS classifiers and the WINGS visual classification (left 
panel) and the automatic classification performed by MORPHOT (right panel). The circle radius is proportional to the number of 
galaxies. 



galaxy masses at high-z, and because it can be used also for 
galaxies with no spectroscopy in the magnitude- "delimited" 

sample. 

We use the equation that lBell fc d e Jong f 200 1') give for 
the Bruzual & Chariot model with a .Salpeter, (,195a) IMF 
(0.1-125 Mq) and solar metallicity, 



logio(M/Lfl) = -0.51 + 1.45 x(B~V) 



(1) 



The total luminosity Lb has been derived from the to- 
tal (S Extractor AUTO) observed B magnitude (|Varela et al.l 
I2009D . corrected for distance modulus and foreground 
Galaxy extinction, and k-corrected using tabulated values 
from (jPoggianti 1997]). The (B — V) color used to calculate 
masses was derived from observed B and V aperture mag- 
nitudes measured within a diameter of 10 kpc around each 

galaxy baricenter, corrected as the total magnitude^ 

Then, we scaled our masses to the more used iKroupal 

l|200ll) IMF adding -0.19 dex to the logarithmic value of the 
meisses. 

Stellar masses for WINGS galaxies observed spectro- 
scopically had been previously determined by fit ting the op- 
tical spectrum (in the range ~ 3600— ~ 7000 A) ([Fritz et al.l 



2010l). with the sp ectro-photometric model fully described in 
Fritz et al] (|2007l ) , and they are in good agreement with the 



masses used in this paper. For a detailed description of the 
determi nation of m a sses a nd a compa rison of different meth - 
ods seelFritz et al.l (|2()10|). F ig. 1 in IVulcani et all (|2010bl ) 
and IValentinuzzi et aLT i 20 id ) . 



2.1.4 Mass-limited sample 

For the mass-limited sample, we rely on spectroscopy to be 
sure that we are using only cluster members. In WINGS, 
photo-z techniques cannot be used to assess cluster mem- 
bership due to the low redshift and to the fact the we have 
at our disposal few photometric bands. Galaxies are consid- 



ered members of a cluster if their spectroscopic redshift lies 
within ±3cr from the cluste r mean redshift, where a is the 
cluster velocity dispersion (|Cava et al.l 12003 ). We use only 
spectroscopically confirmed members of 21 of the 48 clus- 
ters with spectroscopy. The clusters used in this analysis are 
listed in Table [T] This is the subset of clusters that have a 
spectroscopic completeness (the ratio of the number of spec- 
tra yielding a redshift to the total number of galaxies in the 
photometric catalog) larger than 50%. We apply a statisti- 
cal correction to correct for incompleteness, weighting each 
galaxy by the inverse of the ratio of the number of spectra 
yielding a redshift to the total number of gal axies in the 
V-ba nd photometric catalog, in bins of 1 mag (|Cava et al.l 
I2OO9I ). 

In each cluster, we exclude the brightest cluster galaxy 
(BCG), defined as the most luminous galaxy of each cluster, 
that has pecul iar properties and could alter the distributions 
(Fasano et all [2010 '). 

Only galaxies lying within O.6-R200 are considered, be- 
cause this is the largest radius covered in the 21 clusters 
considered (except for A 1644 and A3266 where the cover- 
age extends to ~ O.5-R200). 

To determine the galaxy stellar mass limit of our sam- 
ple, we compute the mass of an object whose observed mag- 
nitude is equal to the faint magnitude limit of the survey, 
and whose color is the reddest color of a galaxy at the high- 
est redshift considered. 

The spectroscopic magnitude limit of the WINGS sur- 
vey is V=20. Considering the distance module of the most 
distant WINGS cluster is ^ 37.5 , and the reddest galaxy has 
a color of {B — V) = 1.2, the magnitude limit corresponds 
to a mass limit Af, — 1O^'*M0, above which the sample is 
unbiased. Adopting this limit, the final sample consists of 
951 early-type galaxies, of which 364 are ellipticals and 587 
are SOs. The corresponding numbers weighted for incom- 
pleteness are: 1469 early-types, 557 ellipticals and 912 SOs. 
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Table 1. List of WINGS clusters analyzed in the mass-limited 
sample, their redshift, velocity dispersion, distance modulus and 
R200- 



Numbers of WINGS galaxies above the EDisCS mass limit 
(M, = 1O^°-^M0, see below) are 594 early-types (920 once 
weighted), 224 ellipticals (341 once weighted) and 370 SOs 
(579 once weighted) (see Table [2]). 



2.1.5 Magnitude- "delimited" sample 

For the magnitude- "delimited" sarnple we use the photo- 
metric data for 76 WINGS clusterslj The clusters used are 
those presented in Table [1] plus those in Table [3] 

We foll o w th e selection criteria proposed by 
iHolden et al.l (|2009l ). They selected a sample of early- 
type gala:xies (ellipticals and SOs) that lie on the red 
sequence (determined with spectroscopic members and 
accepting all gala xies lying within 2a from the sequence 
- for details see iMei et al] |2009| ). At low-z, they used 
only spectroscopically confirmed members. At high-z they 
considered all red sequence galaxies in the photometric 
catalogue except those that are interlopers confirmed by 
spectroscopy. Moreover, they selected galaxies within a 
magnitude range, taking into account passive evolution: 
-19.3 > Mb + 1.208z > -21. Finally, they considered 
only galaxies within 2_R2()o/7r of the cluster center. They 
computed e llipticitics using the results from GALFIT 
l|Peng et al.l |2002) and a dopte d visual morphologies from 
the literature (iDressletl 1980l for the s ample at low-z, 
iDesai et al.ll200'i1 and lPostman et al.ll200^ for the sample at 
high-z) . 

For our WINGS magnitude- "delimited" sample, to be 
strictly consistent with what we do for the EDisCS data- 
set, we consider all galaxies in the photometric catalogue, 
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R200 






{km s~^) 


(mag) 


kpc 


A85 


0.0521 




36.83 


2.54 


A133 


0603 


810±78 


37.16 


1.95 


A147 


0.0447 


UUU^ -LO 


36 48 


1.61 


A160 


0.0438 


'jyj J — Lijvj 


36.44 


1.36 


A168 






36.49 


1.22 


A193 


0.0485 




36.67 


1.83 


A311 


0.0657 


NULL 


37.35 


0.00 


A376 


0.0476 


852di49 


36.62 


2.06 


A548b 


0.0441 




36.45 


2.05 


A602 


0.0621 


720di73 


37.22 


1.73 


A671 


0.0507 




36 77 


2.19 


A780 


n O'^fl^ 


734dil0 


37.01 


1.26 


A1291 


n (T^nQ 

u. uouy 


'i^C'ZIZrty 


36.77 


1.04 


A1668 


0.0634 


649±57 


37.27 


1.56 


A1736 


0.0461 


SKQ-|_fiO 


36.55 


2.06 


A1795 




72c; -LKQ 

( • ' ' '■ 1 


37.27 


1.74 


A1831 


0.0634 




37.27 


1.30 


A1983 


0.0447 




36.48 


1.28 


A1991 


0.0584 




37 08 


1.44 


A2107 


0.0410 




36.29 


1.44 


A2124 


0666 


801it64 


37.38 


1.92 


A2149 


0.0675 




37.41 


0.85 


A2169 


0.0578 


> n in 1 '-HI 


37.06 


1.22 


A2256 


0.0581 


1 273-1-64 


37.07 


3.06 


A2271 


0.0584 


504itl0 


37.08 


1.21 


A2457 


0.0584 




37.08 


1.40 


A 2572a 


0390 


631itl0 


36.18 


1.53 


A2589 


0.0419 


111 1 f 1(1 


36.34 


1.98 


A2593 


0.0417 


701it60 


36.33 


1.70 


A2622 


0.0610 


696-1-55 


37.18 


1.67 




0.0548 




36.94 


1.51 


A2657 


0.0400 


QOI _|_OQ 
vJO-L _l_OtJ 


36.23 


0.92 




0562 


NULL 


37.00 


0.00 


A2717 


0.0498 


■ l< 1. 1 1 ■ 1 j1 


36.73 


1.34 


A2734 


0.0624 


KKK -1-42 


37.23 


1.33 


A3164 


0.0611 


NULL 


37.19 


0.00 


A3497 


u.uuou 


726it47 


37.43 


1.74 


A 3528a 


0535 


<J<Jxy _HJrt 


36.89 


2.17 




\j .yjooo 




36.89 


2.08 


A 3530 


0.0544 


■ M 1. > 1 ■ 1 


36.92 


1.36 


A3532 




K21 -l-^'H 


36.97 


1.50 


A3558 


0.0477 


Qi c;-|-KO 


36.63 


2.21 


A3667 


u.uoou 




36.87 


2.39 


A3716 


0.0448 




36.49 


2.07 


A3880 


0.0570 


763±65 


37.03 


1.84 


A4059 


0.0480 


715±59 


36.64 


1.73 


IIZW108 


0.0483 


513±75 


36.66 


1.24 


MKW3s 


0.0444 


539±37 


36.47 


1.30 


RX0058 


0.0484 


637±97 


36.66 


1.54 


RX1022 


0.0548 


577±49 


36.94 


1.39 


RX1740 


0.0441 


582±65 


36.45 


1.41 


Z1261 


0.0644 


NULL 


37.31 


0.00 


Z2844 


0.0503 


536±53 


36.75 


1.29 


Z8338 


0.0494 


712±60 


36.71 


1.72 


Z8852 


0.0408 


765±63 


36.28 


1.86 


le 3. List of 


additional 


WINGS 


clusters 


used for 



magnitude- "delimited" sample, their redshift, velocity dispersion, 
distance modulus and -R200- 



^ A3562 was excluded due to bad V-band seeing. 
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WINGS 






EDisCS 






M, > 10 




M, > 10 




mag 


M, 3s IO^O-^Mq 


mag 




No6» 


N„ 




N„ 


N 


N 




ellipticals 


364 


557 


224 


341 


580 


145 


101 


SOs 


587 


912 


370 


579 


914 


61 


43 


early-types 


951 


1469 


594 


920 


1494 


206 


144 



Table 2. Number of galaxies in the mass-limited and in the magnitude- "delimited" samples. For the mass-limited sample, for WINGS 
both observed numbers and numbers weighted for spectroscopic incompleteness are given. 



excluding those that are non-members based on the spec- 
troscopy. Contamination on the red sequence at low-z is min- 
imal, and we have checked that the results remain the same 
using only spectroscopic members corrected for complete- 
ness. We exclude from our analysis galaxies located outside 
of -R200, to be consistent with what we do at high-z (see 
§2.2.2). 

We then select only gala x ies ly ing within 2a from the 
red sequence. Like iMei et al] l|2009l '). we define the red se- 
quence using only spectroscopic members and we build a 
color-magnitude diagram for each cluster. To do this, we use 
the observed B and V aperture magnitude measured within 
a diameter of 5 kpc around each galaxy baricenter and the 
total V SExtractor AUTO magnitude, both corrected for 
distance modulus and foreground Gal axy extinc t ion, a nd k- 
corrected using tabulated values from iPoggiantil (|l997h . For 
those clusters for which spectroscopy was not available, we 
define the red sequence using the photometry of morpholog- 
ically selected early-type galaxies. We determine the slope 
and the dispersion of the red sequence in the color- magni- 
tude diagram perfoming a weighted least-square-fit on the 
data, giving less weight to the outliers and reiterating 10 
times to have a better determination of the parameters. 

We include galaxies with -19.3 > Mb + I.2O82 > -21, 
where Mb is the magnitude derived fr om the total (SEx - 
tractor AUTO) observed B magnitude (jVarela et al.ll2009^ ■ 
corrected as the color. In this way we automatically exclude 
the BCGs. Finally, we consider only galaxies that are ellip- 
tical and SO following our morphological classification. 

Our final magnitude-" delimited" sample consists of 580 
ellipticals and 914 SOs, for a total of 1494 early-type galaxies 
(see Table [2|. 



2.2 High-z sample: EDisCS 

The multi- wavelength photometric and spectroscopic sur - 
vey of distant clusters named EDisCS (jWhite et al.l I2OO5I ') 
has been developed to characterize both the clusters them- 
selves and the galaxies within them. It observed 20 fields 
containing galaxy clusters at 0.4 < 2 < 1. 

Clusters were drawn from the Las Campanas Distan t 
Cluster Survey (LCDCS) catalog (|Gonzalez et all bOOlT ). 
They were selected as surface brightness peaks in smoothed 
images taken with a very wide optical filter (~ 4500 — 7500 
A). The 20 EDisCS fields were chosen among the 30 highest 
surface brightness candidates, after confirmation of the pres- 
ence of an apparent cluster and of a possible red sequence 
with VLT 20 min exposures in two filters (jWhite et al.l 
I2OO5I ). 

For all 20 fields, EDisCS has obtained deep optical 



multiband photometry with F0RS2/VLT (|White et all 
120051^ and near-IR p hotom etry with SOFI/NTT 



( Aragon-Salamanca et all I2OI0I ). Photometric redshifts 
were measured u sing both optical and inf rared imaging (see 
IPello et al.ll2009l and lRudnick et al.|[2009l for details). They 
were computed for every object in the EDisCS fields using 
two independent codes , a modified version o f the publicly 
av ailable Hyperz c ode (jBolzonella et al.ll2000l ) and the code 
of Rudni ck et al.l ll200ll) with the modifications presented 
in Rudn ick et al.l (|2003l ). Photo-z membershi p (see also 
[Pe Lucia et al.l 12004 and iDe Lucia et al.l |2007| for details) 
was established us ing a modified ve r sion o f the technique 
first developed in iBrunner fc LubinI (|2000l ). in which the 
probability of a galaxy to be at redshift z {P{z)) is inte- 
grated in a slice around the cluster redshift to give Pdust for 
the two codes. A galaxy was rejected from the membership 
list if Pciust was smaller than a certain probability Pthresh 
for either code. The Pthresh value for each cluster was 
calibrated from EDisCS spectroscopic redshifts and was 
chosen to maximize the efficiency with which spectroscopic 
non-members are rejected while retaining at least ~ 90% 
of the confirmed cluster members, independent of their 
rest-frame (B-V) color or observed (V-I) color. In practice 
it was possible to choose thresholds such that this criterion 
was satisfied while rejecting 45%-70% of spectroscopically 
confirmed non-members. Applied to the entire magnitude 
limited sample, these thresholds reject 75%-93% of all 
galaxies with Itot < 24.9. A posteriori, it was verified that 
in the sample of galaxies with spectroscopic redshift and 
above the mass limit described below, 20% of those galaxies 
that are photo-z cluster members are spectroscopically 
interlopers and, conversely, only 6% of those galaxies 
that are spectroscopic clustermembers are rejected by the 
photo-z technique. 

Deep spectr oscopy with F0RS2 / VLT was obtained for 
18 of the fields (|Hallidav et al.ll2004l : iMilvang -Jensen et al.l 
[2008?). Spectroscopic targets were selected from I-band cat- 
alogs, producing an essentially I-band selected sample with 
no selection bias down to 7 = 22 at 2 ~ 0.4 — 0.6 and 7 = 23 
at z ~ 0.6 — 0.8 (|Hallidav et al. l l2004l : lMilvang -Jensen et al.l 
l2008h . Typically, spectra of more than 100 galaxies per field 
were obtained. 

ACS/HST mosaic imaging in Fil4 W of 10 of the high - 
est redshift clusters was also acquired l|Desai et al.ll2007l l. 
covering with four ACS pointings a 6.5' x 6.5' field with 
an additional deep pointing in the center. This field cov- 
ers the 7?2oo of a ll clusters, except for 1232.5-1250 where it 
reaches 0.57?2oo l|Poggianti et al.ir2006l '). The 7?2oo values for 
ou r structures are compu ted from the velocity dispersions 
bv lPoggianti et al.l (|2008l ). 
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cluster name 



o" -R200 
{kms~^) Mpc 



cl 1040.7- 
cl 1054.4- 
cl 1054.7- 
cl 1103.7- 
cl 1138.2- 
cl 1216.8- 
cl 1232.5- 
cl 1354.2- 



1155 
1146 
1245 
1245 
1133 
1201 
1250 
1230 



0.70 
0.70 
0.75 
0.62 
0.48 
0.79 
0.54 
0.76 



418 
589 
504: 
336 
732 
1018 
1080 
648^ 



-55 

-46 
-78 
-70 
113 
65 
-1-36 
40 
72 
76 
+73 
-77 
-119 
-89 
105 
110 



0.70 
0.99 
0.82 
0.41 
1.41 
1.61 
1.99 
1.08 



Table 4. List of EDisCS clusters analyzed in this paper, 
with cluster name, r edshift, velocity dispersion and -R200 (from 
Hallidav et"ai]|2004l: fMilvang-Jensen et al.ll2008l ; iPoggianti et al.l 
2008t) . 



2.2.1 Morphologies, ellipticity measurements and galaxy 
stellar masses 



Morphologies are discussed in detail in iDesai et al.1 l|2007l ). 
The morphological classification of galaxies is based on the 
visual classification of HST/ACS F814W images sampling 
the rest-frame ~ 4500 - 5500 A range, similarly to WINGS. 

The det e rmina tion of ellipticities is presented in 
ISimard et al.l ('2009) . They have been estimated using the 
tool G IM2D (Galaxy IMa ge 2D) version 3.2, a fitting pro- 
gram ()Simard et al.ll2002l l that performs a detailed surface 
brighteness profile analysis of galaxies in low signal-to-noise 
(S/N) images in a fully automated way. In this paper, we use 
the ellipticities derived fitting every source in the HST/ACS 
images with a single Sersic fit model. Since the shape of 
the PSF on the HST/ACS images varies significantly as a 
function of position, spatially-varying PSF models for the 
EDisCS cluster images were constructed. 

In this analysis, we consider only 8 of the 10 EDisCS 
clusters for which //ST images are available. In fact, elliptic- 
ity measurements are not available for cl 1227.9-1138 and 
no galaxies of cl 1037.9-1243 enter our final samples (see 
below the selection criteria). 

For EDisCS galaxies, we use stellar masses estimated 
using the same relation w e use for the WINGS dataset, 
hence again following the iBell fc de Jon3 ll200ll) method 
and then converting masses to a iKroupal ()200ll ') IMF. To- 
tal absolute mag nitudes are derived from photo-z fitting 
JPello et al.ll2009|). rest-fram e lu minosities have been derived 
using iRudnick et all (l2003h andlRudnick et aP ()2006h meth- 
ods and presented in iRudnick et al.l l|2009l ). Stellar masses 
for spectroscopic members were also estimated using the 
kcorrect tool (iBlanton fc Roweis ,2007. ) that yields masses 
in agreement with those used in this paper. For a detailed 
discussion of our mass esti mates and of the consis tency be- 
tween different methods see IVulcani et all l|2010bh . 



2.2.2 Mass-limited sample 

For the EDisCS mass-limited sample we use all photo- 
z members, following the membership criteria described 
above. 



http: / / cosmo.nyu.edu / mbl44 /kcorrect / 



The choice to use the photo-z membership instead of 
spectroscopically confirmed members is dictated by the fact 
that otherwise the number of galaxies would be low, not 
allowing a statistically meaningful analysis. 

Moreover, the spectroscopic magnitude limit ranges be- 
tween 1=22 and 1=23 depending on redshift, and the corre- 
sponding spectroscop ic stellar mass limit is M = lO^^'^M© 
jVulcani et al.l 120101 ) . The photo-z tecnique allows us to 
push the mass limit to much lower values than the spec- 
troscopy. We adopt a conservative magnitude completeness 
limit for the EDisCS photometry equal to 7 ^ 24 (though 
the completeness remains very high to magnitudes signifi- 
cantly fainter than / = 24, White et al. 2005). We consider 
the most distant cluster, cl 1216.8-1201, that is located at 
2 ~ 0.8 and determine the value of the mass of a galaxy 
with an absolute B magnitude corresponding to / = 24, and 
a rest-frame color (B — V) ~ 0.9, which is the reddest color 
of galaxies in this cluster. In this way, the EDisCS mass com- 
pleteness limit based on photo-z's is M« = 10^^'^ Mq. This 
is the mass limit we a dopt for our analysis. As we discuss in 
I Vulcani et aI](|2010b^ ■ spectroscopic and photo-z techniques 
give very consistent results for the galaxy mass functions 
in the mass range in common. Also comparing the elliptic- 
ity distribution determined using spectroscopic and photo- 
z data down to the spectroscopic mass limit, we find that 
they are not statistically different: a Kolmogorov-Smirnov 
test cannot reject the null hypothesis that the distributions 
are drawn from the same parent distribution with a prob- 
ability of ~ 22% (for details on the K-S test see This 
gives additional support to our choice to use photo-z data. 

As for the WINGS mass-limited sample, both BCGs 
and all galaxies at radii greater than r = 0.6J?200 have been 
excluded from the analysis. Table |4] presents the list of clus- 
ters used and some relevant values. 

The final mass-limited EDisCS sample of galaxies with 
a measured ellipticity for M« ^ IO^'^'^Mq consists of 206 
early-type galaxies, 145 of which are classified as ellipticals 
and 61 as SOs (see Table [2j. 

2.2.3 Magnitude- "delimtted" sample 

For the magnitude - "delimited" sample at high-z, to follow 
the same criteria of lHolden et al.l(|2009l ). we do not consider 
the photo-z membership, but we exclude only those galaxies 
that have been identified spectroscopically as non-members. 
Then, we use all early-type galaxies within 2a of the red 
sequence and with -19.3 > Mb + 1.208z > -21. 

To d e termin e the red sequence of each cluster, as 
iMei et all lj2009l }_ did, we use only spectroscopic members 
of ou r clusters (jHallidav et al ] I2OO4I : iMilvang- Jensen et al.l 
l2008h . We build color-magnitude diagrams using the R — I 
color (that corresponds to ~ S — 1/ in WINGS). Only for cl 
1232.5-1250 we use the V — I color because the R band is 
not available. 

Similarly to what we do for WINGS, we determine the 
red sequence by performing a weighted least-square-fit of 
our data. However, since the resulting red sequences are not 
always reliable, for all clusters but cl 1232.5-1250, for which 
we use a different color, we determine a mean slope using 
only cl 1216.8-1201 and cl 1054. 4-1146, two clusters located 
almost at the same redshift for which the red sequence is well 
defined from the spectroscopy, and adopt this slope for all 
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clusters. The mean dispersion is determined averaging the 
dispersion of all clusters except cl 1354.2-1230 that has been 
excluded because, having too few points, it would give too 
small a value of the dispersion. Then, we determine sepa- 
rately the red sequence of each cluster, using the slope and 
the dispersion just determined and finding the most appro- 
priate value of the intercept. Obviously, the red sequence 
for cll232. 5-1250 is determined separately using its spectro- 
scopic data. 

Subsequently, we consider only galaxies within -R200 (as 
for WINGS), instead of 2i?20o/7r as Holden et al. (2009) did, 
to improve the statistics. 

Our final sample consists of 101 ellipticals and 43 SOs, 
for a total of 144 early-type galaxies (see Table (2)1 . 



3 RESULTS: THE ELLIPTICITY EVOLUTION 
IN MASS-LIMITED SAMPLES 

In this section, we analyze the ellipticities of galaxies in our 
mass limited samples. 



3.1 Ellipticity and SO/E number ratio as a 
function of galaixy stellar mass 

Figure [5] shows the trend of ellipticity as a function of galaxy 
stellar mass for early-type, elliptical and SO galaxies for the 
WINGS and EDisCS samples, above their mass complete- 
ness limits. We compute the median values of ellipticity both 
over the whole mass range (green dashed lines) and in mass 
bins of 0.4 dex (red solid lines). 

At both redshifts, the trend of the ellipticity of all 
early-type galaxies together clearly depends on galaxy mass. 
Less massive galaxies tend to have ellipticities that extend 
to much higher values compared to higher mass galaxies 
which populate only the lower end of the range. Consider- 
ing only ellipticals, the trend is much less striking, though 
still present over the whole mass range in WINGS, while a 
drop is observed only above M* ~ IQ^^-^M© in EDisCS. The 
median ellipticity of SOs shows no clear trend with galaxy 
mass, at least below M. ~ IQ^-^M© in WINGS. At high 
masses, both the apparent fall in WINGS and the rise in 
EDisCS may be simply due to low number statistics. 

We note that ellipticals always have an ellipticity lower 
than 0.6 and mostly below 0.4, while SOs cover a wider 
range of ellipticities, with the majority being concentrated 
at high values of ellipticities, above 0.4. Furt hermore, ellip- 
ticals r each higher mass values than SOs (see lVulcani et abl 
l|2010bD for the mass distribution of ellipticals and SOs in 
these samples, see also §5). 

Clearly, the strong trend of ellipticity with mass ob- 
served in early-types is due both to the trend of ellipticity 
of elliptical galaxies with mass, and, mostly, to the fact the 
ellipticals and SOs are found in different proportions at dif- 
ferent masses: SO galaxies, with their average higher elliptic- 
ities, become more frequent going to lower masses. Figure [6] 
shows the ratio of the number of SO to elliptical galaxies at 
different masses. In WINGS (left panel), the SO/EU ratio 
strongly depends on mass: at higher masses there are pro- 
portionally more elliptical galaxies than at lower masses. 
In the highest mass bin, the ratio drops to ~0, indicat- 
ing that there are almost only elliptical galaxies, while at 




Figure 6. Ratio of the number of SO/Ell galaxies at different 
masses both for WINGS (loft panel) and for EDis CS (right pane l) 
in our mass-limited sample. Errors are binomial llGehrelslll986h . 



M. - 10'"-'^ M0 SOs are twice as numerous as ellipticals. In 
contrast, in EDisCS (right panel) we find that the trend is 
almost flat up to A/, ~ IO^^'^A/q and SOs are less than half 
of the ellipticals. 

EDisCS clusters are seen at an epoch prior to the build 
up of the SO cluster population and Figure |6] clearly shows 
that such build-up occurs mainly at masses below 10" Mq. 



3.2 The evolution of the median ellipticity and of 
the ellipticity distributions 

Table[5]summarizes the median values of ellipticities for both 
samples over the whole range of masses. Errors are estimated 
using the bootstrap resampling method. We adopt these es- 
timates because we want to characterize the errors on the 
medians and not the dispersion of the points around the me- 
dian value (that is the standard deviation). In WINGS, the 
choice of the mass limit (Af, ^ lO^'^Af© or M, ^ 10^°-^ Mq) 
does not alter the final results. 

Comparing low- and high-z, the median ellipticity of SOs 
is compatible within the errors at the two redshifts, while 
it slightly changes with redshift for ellipticals, and more no- 
ticeably for early-types. In particular, it slightly decreases 
going to the current epoch for ellipticals, while it clearly in- 
creases for the early-types. This raise for early types isdue 
to the fact that, as shown in Figure |6] the fraction of SOs 
increases at low-z, mainly in the low mass range. Since SOs 
are more flattened than ellipticals, the median ellipticity of 
early- types shifts to higher values at low redshift. As for the 
evolution of the median of elliptical galaxies, this will be 
discussed later in this section. 

We now compare the high- and low-z ellipticity distri- 
bution, 

to see if it evolves. At both redshifts, we consider 
only galaxies above the common mass limit that is Af* = 



10 Vulcani et al. 



0.8 



0.6 



0.4 - 



-I—' 



.0.2 -: 







_0.8 
^ 0.6 
0.4 



0.2 - 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

early 
■ ' . WINGS. 




I I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



early 
EDisCS. 




I 1 1 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 



I I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I I 

ellipticals 
WINGS 




I I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 I 



ellipticals 
EDisCS 




I I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I I 



I I I 1 1 1 I I 1 1 1 I I 1 1 1 I I 1 1 1 I I 

: ■ SOs ~. 
' .WINGS . 




I I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 I 



SOs 

EDisCS 




I I I 1 1 1 I I 1 1 1 I I 1 1 1 I I 1 1 1 I I 



10 10.5 11 11.5 12 10 10.5 11 11.5 12 10 10.5 11 11.5 12 

Log Mq 

Figure 5. EUipticity vs mass in tiie mass-limited sample. Black points: galaxies. Red solid lines: median and 1 u values, where a is 
the rms, estimated in mass bins. Green dashed lines: median and 1 a values computed over the whole mass range. Top panels: WINGS 
data. The WINGS median are corrected for spectroscopic incompleteness. Left: early-type galaxies, Central: elliptical galaxies. Right: SO 
galaxies. Bottom panels: EDisCS data (panels as for WINGS). 





WINGS 


EDisCS 




M*/M0 ^ lO^-S M,/Mq ^ IQio a 


M^/Mq > 1010-2 


ellipticals 


0.190±0.011 0.179±0.011 


0.220±0.011 


SOs 


0.454±0.014 0.462±0.015 


0.493±0.032 


early-types 


0.336±0.012 0.328±0.016 


0.265±0.013 



Table 5. EUipticity median values for both mass-limited samples with errors defined with bootstrap resampling. For WINGS, medians 
are computed taking into account the weights. For WINGS, also values above the EDisCS mass limit are given. 



We build both the cumulative distributions and his- 
tograms (in bins of ellipticity equal to 0.05) for each class 
of galaxies analyzed. For WINGS, both of them take into 
account the spectroscopic completeness weights. 

Figure [7] shows how the ellipticity distribution of early- 
type galaxies evolves with redshift. As expected given the 
evolution of the median ellipticity, there are proportion- 
ally more galaxies with higher ellipticities at low- than at 
high-z, indicating that low-z early-type galaxies are on av- 
erage more flattened. The overall WINGS distribution (blue 
dashed lines) is quite flat, in particular at intermediate 



values of ellipticity. Instead, the ellipticity distribution of 
EDisCS early-type galaxies (black solid lines) shows a peak 
around e ~ 0.2. 

To quantify the differences between the two distri- 
butions, we perform a Kolmogorov-Smironov (K-S) testQ 

^ The standard K-S, in building the cumulative distribution, as- 
signs to each object a weight equal to 1. Instead, our WINGS data 
are characterized by spectroscopic completeness weights. So, we 
modified the test, to make the relative importance of each galaxy 
in the cumulative distribution depend on its weight, and not being 



The evolution of early-type galaxies in clusters 11 



early-types 




0.4 0.6 
e{elliptlcity) 

Figure 7. Comparison of ellipticity distributions of early-type 
galaxies in the mass-limited samples. Medians and bootstrap er- 
rors are also indicated. Black solid lines: EDisCS; blue dashed 
lines: WINGS. Top panel: cumulative distributions of ellipticity. 
Pk—S is the probability that the two distributions are drawn from 
the same parent distribution. Numbers in brackets are the number 
of galaxies in the considered samples. Bottom panel: histograms 
in bins of 0.05 dex, normalized to 1. 



Throughout this paper, we will consider significantly dif- 
ferent two distributions if the K-S test gives a probability 
< 5%. For early-type galaxies, the K-S test allows us to ex- 
clude the similarity of the ellipticity distributions at the two 
redshifts, giving a probability of ~ 0.04%. 

Next, we investigate whether the observed differences 
are due to an evolution of the ellipticity distribution of el- 
lipticals, of SOs, or both. Figure [8] shows that for ellipticals 
results are ambiguous: the K-S test, giving a probability of 
5.84% of similarity of the distributions, is not strictly con- 
clusive. In contrast, for SOs the distributions are compatible 
with being similar (Pa'-s ^20%). EDisCS SOs are very few 
and this could influence the results of the K-S test; however, 
the cumulative distributions appear to resemble each other, 
indicating that the result of the K-S test should be reliable. 

We wish to go deeper into our analysis, trying to under- 
stand if the K-S results are confirmed and above all if they 
are driven by a different shape of the distributions or simply 
by a different location of the two populations. To do this, 
we perform two other statistical non-parametric tests (i.e. 
they do not assume the normal distribution) which make no 
assumptions about the distributions of the populations. 

In Appe ndix E|we present the d e tailed analysis of the 

iMosed (|l963l ) and lMann fc WhitnevI l|l94'<1 ) tests. In sum- 
mary, although the K-S test is inconclusive, from these ad- 



ditional tests it emerges that some mild differences exist in 
the ellipticity distribution of ellipticals at high- and low-z. 



3.2.1 Round ellipticals at low-z 

Inspecting the histograms in Figure |8l it is clear that, at 
both redshifts, there are no ellipticals with e J5 0.6 and both 
distributions are peaked around e ~ 0.2 — 0.3. It seems that 
the greatest differences are confined to the extremes of the 
distributions: in the highest ellipticity bins there are pro- 
portionally more EDisCS elliptical galaxies than WINGS', 
and more noticeably, in the first bin there are proportionally 
more WINGS galaxies with e ^ 0.05 than EDisCS galaxies. 
It is important to stress that this goes in the opposite direc- 
tion of what could be expected from morphological classifi- 
cation biases at high-z: face-on SOs would be systematically 
mistaken for ellipticals more frequently at high- than at low- 
z. Analyzing more accurately the first bin (plot not shown), 
we find that it is mainly dominated by a second peak around 
e ~ 0.03 that instead is not detected in the EDisCS distribu- 
tion. At the moment, we are not able to explain why at low 
redshift there is an exceeding population of rounder ellipti- 
cals compared to high-z; anyway we think it is real, having 
accurately checked both morphologies and ellipticities for 
those galaxies 

We have also visually inspected the ellipticity profiles 
of these round WINGS ellipticals. For a few of them, the 
profile is altered by crowding, for some other the twist of the 
isophotes is very marked and the ellipticity value strongly 
changes with radius. However, 60% of the analyzed galaxies 
really have a very low ellipticity at all radii. Analyzing and 
understanding this population goes beyond the scope of this 
work and it will be discussed in a forthcoming paper. Here, 
we just notice that this populatio n exists and that it could 
be the result of d ry merger events (|van Dokkum et"al ] | 19991 : 
iTran et al.ll2005l ') which might preferentially result in round 
galaxies building over time. 

If elliptical galaxies with e < 0.05 are excluded, the 
WINGS and EDisCS ellipticity distributions for ellipticals 
are indistinguishable (the K-S test gives a probability > 20% 
that the poopulations are drawn from the same parent dis- 
tribution). 



3.3 What drives the evolution of the ellipticity 
distribution of early-type galaxies: the 
evolution of the galaxy mass distributions, or 
the evolution of the relative proportions of 
ellipticals and SOs? 

To summarize the most important points, in the previous 
sections we have found that, for mass-limited samples, the el- 
lipticity distribution of early-type galaxies strongly changes 
with redshift. The ellipticity distributions of ellipticals and 
SOs on the whole do not evolve significantly, although for 
ellipticals there is a non negligible shift of the medians of 



fixed to 1. In the following, we will always use this modified K-S 
test. Obviously, when using photo-z's all galaxies have a weight 
equal to 1, and using the modified test is equivalent to using the 
standard one. 



^ We find a similar peak also in our low-z magnitude- "delimited" 
sample (see §4) a nd also in the low-z sample analyzed by 
iHolden et ahl ||2009D (see §6), so it can neither due to the adopted 
selection criteria nor to a bias in our samples. 
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Figure 8. Comparison of ellipticity distributions of elliptical (left panel) and SO (right panel) galaxies in the mass-limited samples. 
The medians and their bootstrap errors are also indicated. Black solid lines: EDisCS; red (in the left panel) and green (in the right 
panel) dashed lines: WINGS. Top panels: cumulative distributions of ellipticity. Pk~s is the probability that the two distributions are 
drawn from the same parent distribution. Numbers in the brackets are the number of galaxies the considered samples. Bottom panels: 
histograms in bins of 0.05 dex, normalized to 1. 



the distributions due to a relative excess of round ellipticals 
at low-z compared to high-z. 

We have found in Figure [5]that above all for early- types 
there is a trend of ellipticity with mass, and that these trends 
are different in EDisCS and in WINGS. M oreover, as we 
have d iscussed at length in a previous paper jVulcani et al.l 
bOlObl ). WINGS and EDisCS have different galaxy stellar 
mass distributions. We have found that the mass distribu- 
tion of each morphological type evolves with redshift and 
that all types have proportionally more massive galaxies at 
high- than at low-z. As a consequence, if we want to under- 
stand the origin of the ellipticity distribution of galaxies of 
different morphological types and at different redshifts, we 
have to try to disentangle the influence of the evolution of 
the mass distribution from the effects of the morphological 
evolution. 

3. 3. 1 The evolution of the galaxy mass distribution 

First, we analyze separately ellipticals and SOs. We perform 
1000 Monte Carlo simulations extracting randomly from the 
WINGS sample a subsample with the same mass distribu- 
tion as the EDisCS sample, separately for the two different 
morphological classes. Each time we extract from WINGS 
the same number of galaxies that are in the EDisCS sample 
(i.e. 145 ellipticals and 61 SOs). For each simulation we de- 
termine the ellipticity distribution, the median value of the 
ellipticity and we perform a K-S test to compare the result 
of the simulation with the EDisCS data-set. Then, taking 
into account all the simulations, we determine both the me- 
dian value of the K-S and the fraction of simulations that 
give a conclusive statement (see Table [6]). 

Figure [9] shows our results for the two morphological 



classes. The plotted WINGS histogram and cumulative dis- 
tribution (dashed lines) are the average of the 1000 Monte 
Carlo simulations. The K-S probabilities and WINGS me- 
dians reported in the plots are the medians of the values 
of all the simulations. Adopting the same mass distribu- 
tion, for elliptical galaxies, the match of the mass distribu- 
tions at high- and low-z makes the ellipticity distribution of 
WINGS galaxies (red dashed lines in the left panel of Fig- 
ure [9| become significantly different from the EDisCS one. 
The K-S test gives conclusive results {Pk-s ^ 5%) in 80% 
of the simulations, with a median value of 1.3%, rejecting 
the null hypothesis of similarity of the two distributions. 
This seems to imply that WINGS ellipticals tend to be on 
average rounder than EDisCS ellipticals of the same mass. 
This effect is mitigated in the observed WINGS vs EDisCS 
distributions (Figure[8]) (yielding a non-significant K-S test) 
by the fact that a) there are proportionally more less mas- 
sive galaxies at low- than at high- z (for details on the mass 
functions see I Vulcani et al" I l2010bl ): and b) low -mass ellipti- 
cals are more flattened (on average) than high-mass ones. 
Hence, the increase in the number of low mass ellipticals 
at low-z largely compensates the existence of a significant 
evolution in the ellipticity distribution of ellipticals with the 
same mass distribution (Figure [9|, and produces an ambigu- 
ous or at best weak evidence for evolution in our analysis of 
§3.2. 

If we neglect galaxies with e < 0.050 for which the most 
oustanding differences are detected (see §3.2), adopting the 



° In this way we wish to compare the whole general distribution 
without being too much influenced by galaxies located in only 
one bin. 
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Figure 9. Ellipticity distribution of elliptical (left panel) and SO (right panel) galaxies assuming for WINGS the same mass distribution as 
EDisCS (cf. with Figure|8]l. Panels and symbols are the same as in Figure[8] The plotted WINGS histogram and cumulative distributions 
are the average of the 1000 Monte Carlo simulations. The Pk—S the median K-S value. 



Pk-S 

< 1% 1% - 5% 5% - 10% > 10% median 



mass-matched simulations 

eUipticals 50% 30% 9% 11% 1.3% 

SOs 1% 6% 4% 89% >20% 

early-types 18% 30% 20% 32% 5.2% 



morphology-matched simulations 
early-types 0% 3% 8% 89% >20% 



Table 6. Results of the Kolmogorov-Smirnov test performed on the 1000 mass-matched and morphology-matched simulations (see text 
in §3.3). Pk—S is tli^ probability that the WINGS and EDisCS distributions are drawrn from the same parent distribution. 



same mass distribution, the median value of the K-S test 
is ~ 8%, indicating that the major differences between the 
ellipticals at high- and low-z is indeed the enhanced pop- 
ulation of round ellipticals at low-z. Turning to SOs (green 
dashed lines in right panel of Fig. [9|, the mass matched 
ellipticity distributions of WINGS and EDisCS remain sta- 
tistically indistiguishable: the K-S test cannot reject the null 
hypothesis, giving a probability ^ 5% in 93% of the simula- 
tions and a median value > 20%. 

Now we wish to test if the different mass distribution at 
different redshifts alters the ellipticity distribution of early- 
types, so we put together WINGS ellipticals and SOs, being 
sure to extract randomly galaxies in order to have the same 
EDisCS mass distribution and maintaining the WINGS mor- 
phological fractions (i.e. ~ 40% ellipticals and ~ 60% SOs 
- see Table [7]). In this way we test whether the observed 
evolution of the ellipticity distribution of early-type galax- 
ies can be entirely explained by the evolution of the mass 
distributions. 

In Figure [To] the K-S gives a probability ^ 5% that the 
two distributions are driven from the same parent distri- 



bution in 48% of the simulations (Table [SJ, and a median 
probability for all simulations of 5.2%. Excluding galax- 
ies with e < 0.05, the K-S results remain similar (median 
Pk-s ~ 4%, Pk-s < 5% in 55% of the simulations, plot 
not shown) . These values suggest that the different mass dis- 
tribution at the different redshifts influences at some level 
the evolution of the ellipticity distribution, even if proba- 
bly it is not the main factor as it cannot fully explain the 
observed evolution. 

3.3.2 The evolution of the morphological fractions 

We now wish to assess the role played by the evolution 
with redshift of the relative morphological fractions. In 
Table [7] we show how much the morphological fractions 
change with time, both for the mass-limited sample and for 
the magnitude- "delimited" one (see We observe that 
while at high redshift ellipticals are more common than SOs 
(~ 70% and ~ 30% respectively), in the Local Universe SOs 
dominate, representing ~ 62% of the early-types. 

To analyze the importance of this evolution, we now 
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Table 7. Relative morphological fractions of galaxies in both mass and mag-(de)limited samples. For the WINGS mass-limited sample, 
also numbers above the E DisCS mass lim it are given. In both cases, both observed and completeness-weighted numbers are listed. Errors 
are binomial, as defined in iGehrelj l ll986l) . 
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Figure 10. EUipticity distribution of early-type galaxies assum- 
ing for WINGS the same mass distribution as EDisCS galaxies 
and maintaing the WINGS morphological mix (cf. with Figure[7ll. 
Panels and symbols are the same as in Figure [7] The plotted 
WINGS histogram and cumulative distributions are the average 
of the 1000 Monte Carlo simulations. The Pk—S is the median 
K-S value. 



ellipticals and SOs) with e < 0.05, since their presence likely 
alters the final results. In fact, increasing the number of 
elliptical galaxies in WINGS from being -40% to ~70% 
of the whole population, the contribution of galaxies with 
e ~ 0.03 is hugely magnified and it would strongly infiuence 
the whole population. For galaxies with e ^ 0.05, Figure [TT] 
shows that, if at both redshifts we had the same fractions 
of ellipticals and SOs, the ellipticity distributions for early- 
type galaxies would be indistinguishable. The K-S test is 
conclusive {Pk-s ^ 5%) only in 3% of the simulations, while 
the median value of the probability is > 20%. 

Importantly, doing the same in the observed distribu- 
tions and excluding all galaxies with e < 0.05 from Fig. 7, 
the low- and high-z early-type ellipticity distributions re- 
main significantly different, with a K-S test probability to 
be drawn from the same parent distribution of only 0.002%. 
Except for the excess of round ellipticals at low-z, the evolu- 
tion of the ellipticity distribution of early-type galaxies can 
be fully explained by the morphological evolution. 

From this whole section we conclude that it is mainly 
the relative contribution of each morphological type to the 
total that is responsible for the evolution of the elliptic- 
ity distribution of early- type galaxies, even if the role of 
the evolution of the mass distribution with redshift is non- 
negligible. Morphology appears to be the most decisive fac- 
tor in the evolution of the ellipticity distribution of early- 
type galaxies: WINGS and EDisCS have a different morpho- 
logical mix and their ellipticity distribution is regulated by 
the different relative proportions of ellipticals and SOs. 



perform a second test, extracting randomly from the 
WINGS data-set a subsample of galaxies with approxi- 
mately the same relative fraction of SOs and ellipticals as 
EDisCS (i.e. in each simulation we extract 70 ellipticals and 
30 SOs from the WINGS sample) and paying attention to 
maintain the WINGS mass distribution. In this way we wish 
to test whether the morphological evolution can account for 
the ellipticity evolution of early-type galaxies, letting the 
mass distribution to naturally evolve. We perform 1000 such 
simulations. 

We then compare the ellipticity distribution of the 
"modified" WINGS sample, to the real EDisCS one. Con- 
sidering all galaxies (also those with e < 0.05), (plot not 
shown) the K-S test is conclusive {Pk-s ^ 5%) in 69% of 
the simulations, while the median value of the probability is 
0.2%, indicating that even assuming the same morphologi- 
cal fraction in the two samples, differences between the two 
ellipticities distributions are still detected. 

Next, we exclude from our analysis those galaxies (both 



4 RESULTS: THE ELLIPTICITY 
DISTRIBUTION IN THE 
MAGNITUDE- "DELIMITED" SAMPLE 

From EDisCS and WINGS, we have also selected a 
magnitude- "deli mited" sample o f gala xies following the cri- 
teria adopted by iHolden et al.l (|2009l ) . as described in §2, 
with the aim to directly compare our finding with their re- 
sults. 

First of all, we wish to check if the different selection 
criteria implicate a change in our findings compared to the 
mass-limited sample. 

In the magnitude- "delimited" sample, we qualitatively 
find the same ellipticity-mass relation we found in Figure [5] 
(plots not shown) for the mass-limited sample: no trend of 
ellipticity with mass for SOs, slight trend for ellipticals, and 
a striking trend for early- type galaxies. 
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Table 8. EUipticity median values for the magnitude- "delimited" samples with errors estimated with bootstrap resampling. iHolden et aTl 
(|2009l'l is the sample we use as comparison for the magnitude- "delimited" sample (for details, see text, 
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Figure 11. EUipticity distribution of EDisCS and WINGS early- 
type galaxies, assuming for WINGS the same fraction of ellipticals 
and SOs of EDisCS and maintaing the WINGS mass distribution 
(cf. with Figure [7] and Figure [TOt . Panels and symbols are the 
same as in Figure [7] Only galaxies with e ^ 0.05 are considered. 



Figure 12. EUipticity distribution of early-type galaxies in the 
magnitude- "delimited" samples. Panels and symbols are the same 
as m Figure [T] Cf . with the mass-limited sample shown in Fig- 
ure [T] 

bility of ~ 12%. However, we note that the median ellipticity 
of SOs is significantly lower at low-zQ 



Values of the median ellipticities of the different mor- 
phological types can be found in Table [8] 

Analyzing the ellipticity distributions in the magnitude- 
"delimited" sample, for early- types (Figure [12} the K-S test 
cannot detect a significant evolution with redshift, as it gives 
a probability of ~ 7.7% that the two populations are drawn 
from the same parent distribution. However, the median el- 
lipticities are different: the median is higher at low-z (see 
also Table 9). 

Figure [TJ] shows separately the ellipticity distribution 
for elliptical and SO galaxies. We find that the elhptic- 
ity of elliptical galaxies evolves noticeably with redshift: in 
WINGS (red dashed lines in the left panel of Figure [TS]) 
there are proportionally more galaxies with low values of 
ellipticity, (e < 0.15), while in EDisCS (black solid lines) 
there is a noticeable peak around e ~ 0.2. The K-S test 
excludes that elliptical galaxies have a common ellipticity 
distribution at the different redshifts, giving a probability 
of ~ 0.41% that they are drawn from the same parent dis- 
tribution. In contrast, the same test cannot distinguish any 
differences between the distributions of SOs, giving a proba- 



4.1 How can results from different samples be 
reconciled? 

Our mass- and magnitude- (de)limited samples give differ- 
ent results. In Appendix[B]we discuss the reasons for the ob- 
served discrepancies, comparing directly the ellipticity dis- 
tributions of the mass- limited and magnitude- "delimited" 
samples for the same type of galaxies at the same redshift. 
Briefiy, the origin of the observed differences in the distri- 
butions lies in the fact that galaxies in the two samples 
are characterized by different properties; in particular the 

^ Performing the Moses test, we can conclude that the early-types 
are unlikely to have the same scale parameter (with a probabil- 
ity of < 5% in 53% of the simulations), while for both ellipticals 
and SOs we cannot exclude a compatibility of the scale param- 
eter, with a probability < 5% in 39% and in 10% of the sim- 
ulations respectively. From the Mann- Whitney test, we cannot 
exclude compatible medians for early-types (with a probability of 
~ 14%), while we can do it for ellipticals and SOs (with a prob- 
ability of 0.07% and 3.4% respectively) (for details on the tests 
see Appendix IAIi . 
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Figure 13. EUipticity distribution of elliptical and SO galaxies in our magnitude- "delimited" samples. Panels and symbols are the same 
as m Figure m Cf. with the mass-limited sample shown in Fig. 181 



magnitude- "delimited" samples are biased and so they are 
not representative of the overall population. In fact, selecting 
early-type galaxies on the red sequence only in the magni- 
tude range —19.3 > Mb + 1.208z > —21, we are loosing (the 
few) most massive galaxies, and a large fraction of the less 
massive galaxies. 



5 COMPARISON WITH LITERATURE 
RESULTS 

We now show that, even adopting the same selection cri- 
ter ia, our results are n ot in agreement with those reached 
bv iHolden et Zl |200i), who investigated the evolution in 
the ellipticity distribution using two magnitude- "delimited" 
samples of cluster early-type galaxies in two redshift ranges. 
Their sample in the local Universe (z = 0.02 — 0.05) con- 
sisted of 10 clusters (for a total of 210 galaxies), while the 
sample in the distant Universe [z — 0.33 — 1.26) consisted 
of 17 clusters (for a total of 487 galaxies) of galaxies with 
HST images. 

With this selection, their main conclusion is that there 
is no evolution neither in the median ellipticity nor in the 
shape of the ellipticity distribution of cluster early-type 
galaxies with redshift from z > 1 to 2 ~ 0. Their median 
ellipticity at 2: > 0.3 is statistically identical with that at 
z < 0.05 and the shapes of the distributions broadly agree. 
Moreover, they find a statistically significant evolution in the 
SOs ellipticity distribution, while they do not detect evolu- 
tion for the ellipticals. 

In Table [S] we compare our median ellipticities with 
theirs. 

Summarizing, comparing the two magnitude- 
"delimited" samples selected in the same way, even if 
we both do not detect an overall evolution in the ellipticity 
distribution for early-types galaxies, we do find an evolution 



in the median which they do not. Moreover, we find a 
significant evolution for ellipticals and no evolution for SOs 
(though again we do find an evolution in the median), while 
they find no evolution for ellipticals and strong evolution 
for SOs. This leads us to believe that the agreement for 
the early-types is not real, but simply due to a particular 
combination of the distributions of ellipticals and SOs. 

We note that both at low and at high redshift, there 
are some clusters that are in common between our and 
their samples: at low redshift there are 4: A 119, A 168, 
A957x, A1983, at high redshift there are 5: cl 1040.7-1155, cl 
1054.4-1146, cl 1054.7-1245, cl 1216.8-1201 and cl 1232.5- 
1250. For these 5 EDisCS clusters in common at high-z, also 



the morphological classifications (|Desai et al.ll2007f l and the 
images used to measure ellipticities are the same. 

So we can use this information to try to understand the 
reason for the discrepancy of the results. 



5.1 Origin of the discrepancies 

We proceed by compar ing separately our WINGS and 
EDisCS samples with the lHolden et al" (|2009 ^) samples, try- 
ing to identify the reason of the observed differences. 



5.1.1 WINGS and Holden et al. (2009) 

In Figur ell4 | we compare our results at 0.04 < z < 0.07 with 
those of lHolden et al.l (|2009l 'l at z < 0.05. While the ellip- 
ticity distributions of early-type galaxies and ellipticals are 
consistent with being similar in the two samples (in both 
cases the K-S test cannot reject the null hypothesis of simi- 
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Figure 14. Comparison between our results and those of iHolden et al at low-z. Top panels: cumulative distributions; bottom 

panels: histograms normalized to 1. Pk—S the probability of the two distr i bution s are drawn from the same parent distribution. 
Left panel: comparison between WINGS (black so lid lines) and Holden et al.l l |2009l ) (blue dashed lines) early-types. Central panel: 
comparison between WINGS (blac k solid lines) andlHoIden etal] (2009?) (red dashed lines) ellipticals. Right panel: comparison between 
WINGS (black solid hues) and iHolden et al.l l|2009^ (green dashed lines) SOs. 



larity of the distributions, giving a probability of 
those of SOs is remarkably different: WINGS SOs (black solid 
lines in the right panel of Figur e [T4)) peak around e ^ 0.5, 
while SOs of lHolden et al.l ()2009[ ) (green dashed lines in the 
same plot) pe ak around e ~ 0.3-0.35, indicating that in the 
iHolden et al.l (2009) sample galaxies are on average rounder 
than in WINGS. The K-S test finds that the distributions 
can be drawn from the same parent distribution with a prob- 
ability of only ~ 2.3%; 

Since we adopted IHolden et al.l (|2009l ) selection criteria, 
if there are differences, they can be due either to differences 
in the morphological classification, differences in the mea- 
surement of the ellipticities, or to variations of some other 
galaxy properties. Here we try to analyze each one of these 
factors. 

Since IHolden et all (|2009l ) draw morphologies from 
iDressled (|l980l ) and in their sample at low redshift there are 
some clus ters in common with WINGS, we select from the 
catalog of lDressleil (ll980D galaxies in common with WINGS 
and see if there are some noticeable differen ces among them. 
On the whole, there are 18 clusters in the |Dressleir(ll980l ) 
catalog that belong also to WING^. 

First of all, we select those galaxies that are early- types 
according to the lDresslej (jTosO) classification and we assign 
to them our measurements of ellipticity. In this way, we can 
compare directly our values of ellipticity with those calcu- 
lated by iHolden et all (|2009l ). Figure [TSl shows the elliptic- 
ity distribution of the two datasets, for early-type, elliptical 
and SO galaxies: the K-S probability is always inconclusive 
{Pk-s > 20%, ~ 8.5%, ^ 20 %) respectively) This is consis- 
tent with the assumption that IHolden et ah! (|2009l ) estimates 
of galaxy ellipticity are compatible with ours. 

Second, we check that using only WINGS clusters classi- 
fied also by Dressier is equivalent to using the whole WINGS 



* This could be a problem linked to th e poor statistic: if we dou- 
ble the number of lHolden et al. I l l2009t) galaxies the test becomes 
conclusive. 

^ For this com parison we include in the analysis al so those clus- 
ters f rom iDres sler ( 1980) that do not enter the iHolden et al.l 
(|2009|) sample, in order to improve the statistics. 



dataset: in fact the elli pticity di s tribut ions for WINGS clus- 
ters in common with bresslej (jigSO") (18 clusters) are m 
agreement with that of our whole sample (76 clusters) (plots 
not shown, K-S probabilities >> 20% for all the morpholog- 
ical types). The 3 distributions (early-types, ellipticals and 
SOs) are very s imilar, indicatin g that the WINGS clusters in 
common with Dressier (19S&) are not a biassed subsample 
of the whole WINGS dataset. 

Since ellipticity measurements do not seem to be re- 
sponsible for the differences between the two datasets, we 
focus our attention on the other possible sources of differ- 
ences, starting with the morphological classifications. 

We find that -27% (79/29 6) of the galaxi es have been 
classified differently from us and lDressleil|l980t ). This corre- 
sponds to the typical agreement between independent clas- 
sifiers, see also i\'2.1\ However, we have checked that they are 
too few to infiuence the overall ellipticity distribution. Even 
re-classifing them and moving them to the other morpho- 
logical class, they do not alter the ellipticity distribution of 
the class in which they have been inserted. 

Moreover, comparing the ellipticity distribution of 
galaxie s belonging to the same morphological class for us 
and for lDresslej (|l980l ). once again we find no significant dif- 
ferences (the K-S test is always largely inconclusive) (plots 
not shown). 

So, the inconsistency is not even linked to the different 
morphological classification, and we have to focus on possi- 
ble biases due to other factors. 

Third, we inve stigate the magnitude distributions of 
iHolden et al] l|2009l) and WINGS samples, to be sure that 
all samples are equally deep. In Figure [16] we compare the 
magnitude distribution of the analyzed samples and we find 
that the magn i tude distribution of the subsample of galaxies 
with iDresslerl lll980l) morphologies (therefore those used by 
iHolden et aiTl2009l ) (red filled histo gram) is very different 
from that of our galaxies, both if we consider only the clus- 
ters that are in common (green filled histogram) and if we 
consider the whole WINGS sample (blue filled histogram). 
Performing a K-S test on the magnitude distributions, we 
find that the whole WINGS sample and the subsample of 
galaxies with IDresslerl (|l980l ) morphologies are drawn from 
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Figure 15. Comparison between tl ie distribution of W INGS cUipticities for galaxies tliat are early-types for IPressleil l|l98Cl ) (black 
solid lines) and the low-z sample of iHoIden et alj ||2009| V Top panels: cumulative distributions; bottom panels: histograms normalized 
to 1. Pk-S is the probab i lity of the two distributions are drawn from the same parent distribution. Left panel: comparison between 
WINGS and iHolden et ahl ||2009| ) (blue dashed lines) early-types. Central p anel: comparis on between WINGS and lHolden et al.1 ||2009| ') 
(red dashed lines) ellipticals. Right panel: comparison between WINGS and lHolden et al.l l|2009l l ( green dashed lines) SOs. 



different parent distributions {Pk~s ~ 2%). Since we are 
following exactly the same selection criteria, the magnitude 
distribut i ons s hould have been similar. It seems th at the 
iDresslej (|l98(]l ). and hence the iHolden et af] (|2009l ). low- 
z sample misses galaxies at fainter magnitudes. Probably 
they are too fa int to have been morphologically classified by 
iDressleJ l|l980l ) who used photographic plates for the classi- 
fication. 

To test whether this bias considerably alters the ellip- 
ticity distribution, we create magnitude matched samples, 
selecting from our WINGS sample a subsample of galax- 
ies with the same magnitude (within ±0.05 mag), the same 
morph ology an d in the same cluster as the sample of lPressleil 
(|l980l ) used bv IHolden eraP (|2009l ). As shown in Figure [HI 
comparing the WINGS magnitude matched simulated sam- 
ple with the Dressier (hence Holden et al. 2009) sample the 
ellipticity distributions are compatible, both for early-types, 
for ellipticals and for SOs (the K-S test can not reject the 
null hypothesis of common origin of the distribution, giving 
a probability always > 20%). 

We conclude that the differences observed between the 
WINGS and the Holden samples in the Local Universe 
are due to the fact that the latter includes on ly those 
galaxi es that were morphologically classified by iDresslei] 
l|l980 ') and does not correspond to a complete sample within 
the adopted magnitud e limits. Therefore, w hen comparing 
WINGS and the low-z IHolden et all (|2009l ) sample we are 
not comparing samples with the same properties, i.e. with 
the same magnitude distribution. 

5.1.2 EDisCS and Holden et al. (2009) 

Now we wish to check if there are some diff erences also be- 
tween EDisCS and the iHolden et af] (|2009l ) sample. 

Having the ellipticities of both samples, we can com- 
pare directly the ellipticity distributions of the two samples 
at high redshift. Figure [18] shows the comparison between 
the EDisCS and Holden high-z ellipticity distributions for 
the different morphological types. We find that there are no 
significant differences {Pk-s ~ 11%, ^ 20% and ^ 20% for 
early-types, ellipticals and SOs respectively). Also compar- 



ing the ellipticity distribution of galaxies belonging only to 
the clusters in common (plots not shown), we can state that 
there are no discrepancies between the samples, as the K-S 
is always largely inconclusive {Pk-s always >> 20%). 

5.2 Conclusions 

To conclude, the d ifferent result s drawn analyzing our sam- 
ples and that of IHolden et al.l l|2009l ) mainly arise from 
the fact that, despite the fact that in principle galaxies 
are selected following the same criteria, actually at low 
redshift they have a different magnitude distribution. The 
IHolden et al.l l|2009t ) low-z sample suffers from incomplete- 
ness at faint magnitudes, likely due to the lack of availability 
of Dressier (1980) morphologies at faint magnitudes. No dif- 
ferences have been detected instead at high redshift between 
Holden's sample and ours. 



6 SUMMARY AND CONCLUSIONS 

In this paper we have analyzed the ellipticity distribution of 
early-types galaxies, and of ellipticals and SOs separately, in 
clusters at z = 0.04 — 0.07 and z ^ 0.4- 0.8. We have taken 
into account both a mass-limited sample and a magnitude- 
"delimited" sample of galaxies. 

• In our mass-limited samples, above the common mass 
hmit (M* ^ lO^'^'^M©) the ellipticity distribution of early- 
type galaxies strongly varies with redshift. This is due to a 
change both of the median and of the shape of the distribu- 
tions with redshift. For ellipticals, no statistically significant 
differences are observed in the high- and low-z distribution, 
even if an evolution of the medians is detected and we ob- 
serve an excess population of round ellipticals at low-z com- 
pared to high-z. Finally, no evolution is observed for SOs. 
The evolution of early-type galaxies is not simply related to 
the different mass distributions at high- and low-z. In fact, 
removing the influence of the mass, the results remain incon- 
sistent. Instead, it is mainly related to the evolution of the 
morphological mix with redshift and hence to the relative 
contribution of ellipticals and SOs at the two epochs. 
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Figure 16. Magnitude distribution of galaxies in different samples. Blue histogram: early-type galaxies in the w hole WINGS sa mple 
(1494 galaxies). Green histogram: early-type galaxies in the WINGS subsamplc of clusters that are in common with lDressleil jlQSOh (547 
galaxies). Red histogram: galaxies early-types in the WINGS sample that are in common with lDressQ lll980l ) (296 galaxies). 



early-types ellipticals SOs 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 



6(ellipticities) e(ellipticities) e(ellipticities) 

Figure 17. EUipticity distribution of magnitude-matched samples (see text for details) for early-types (left panel), ell ipticals ( c entra l 
panel) and SOs (right panel). Panels and symbols are as usual. Black solid lines represent WINGS galaxies present in iDresslerl lll98Cll ) 
catalog; coloured dashed lines represent WINGS mock distributions. 



• As mentioned in the previous point, in our low-z sam- 
ple, we find a population of very round (e ^ 0.05) elliptical 
galaxies that is less conspicuous at high-z. This population 
seems real and not due to selection effects or measurement 
problems. 

• In our magnitude- "delimited" sample, for early-types 
and SOs the evolution is not evident (though the medians of 
both early types and SOs change with z) , while for ellipticals 
we have found a change of the distribution with redshift. 

• The observed differences between the mass-limited sam- 
ple and the magnitude- "delimited" one can be due to the 
different mass distribution of the two samples: in fact in the 
magnitude- "delimited" samples we are loosing some galaxies 
that enter the mass-limited one, both at high but especially 
and more importantly at low masses. 

• Our magnitud e-" delimited" result s are not in agree- 
ment with those of iHolden et al.l (|2009l ). who also analyzed 
a magnitude- "delimited" sample of early- types belonging to 



the red sequence. The main reason of the observed discrep- 
ancy is that, despite galaxies being selected following in 
principle the same criteria, in practice the two low-z sam- 
ples have a d iffere nt magnitude distribution because the 
iHolden et al.l (|2009| ) sample suffers from incompleteness at 
faint magnitudes. 
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Figure 18. Comparison between our results and those of iHolden et al at high-z. All clusters are used in both samples. Top 

panels: cumulative distributions; bottom panels: histograms normalized to 1. Pk-S is the probability of the two distribut ions are drawn 
from the same parent distribution. Left panel: comparison between EDisCS (black solid lines) and iHolden et al.l | |200S|) (blue dashed 
lines) early- types. Central panel: comparison between EDisCS (black solid lines) and IHolden et al. I l l200g^ (red dashed lines) ellipticals. 
Right panel: comparison between EDisCS (black solid lines) and lHolden et al.l l|2009l ) (green dashed lines) SOs. 



APPENDIX A: ADDITIONAL STATISTICAL 
TESTS 

In § 13.21 we characterize the evolution of the ellipticity dis- 
tribution performing the K-S test. Anyway, since this test is 
as general as possible, we wish to go deeper into our anal- 
ysis, trying to understand if the K-S results are confirmed 
and above all if they are driven by a different shape of the 
distributions or simply by a different location of the two 
populations. To do this, we perform two other statistical 
non-parametric tests (i.e. they do not assume the normal 
distribution) which make no assumptions about the distri- 
butions of the populations. 

We use thelMosej (1196.1) test to check the equality of the 
scale parameter, taking into account that each population 
has a difi'erent median. This test is very useful to compare 
the shape of two distributions and to evaluate their disper- 
sion. Following this procedure, we subdivide each population 
into a certain number of groups, each one containing 10 ob- 
servations. For each group, we compute its average and the 
sum of the residuals. Then we put all together the residu- 
als of the two populations, paying no attention which one 
each value belongs to, and we sort them. Afterwards, we 
sum the rank of each population separately and we compare 
the sums. If they are very different, the probability that the 
populations are drawn from the same parent distribution is 
very small. Since this test requires to consider randomly a 
subsample of the observations of the populations, we repeat 
the test 1000 timesl3 It emerges that early-type galaxies at 
different redshifts are unlikely to have the same scale param- 
eter (with a probability < 5% in 85% of the simulations), 
while both ellipticals and SOs show a high compatibility of it 
(with a probability < 5% only in 14% of the simulations in 
both cases), suggesting that the shapes of the distributions 
are similar at different redshifts. 

Then, to test if there could be a shift in location be- 
tw een the populations, we adopt the U-statistic proposed 
bv lMann fc Whitney! l| 19471) . This allows us to assess if there 

Since this test is based on random samples, we do not take 
into account the WINGS's weights. 



are differences in the median values, regardless of the choice 
of the errors adopted to characterize the medians. This pro- 
cedure requires to rank all the values, without regard to 
which population each value belongs too Similarly to what 
we did for the Moses test, we sum the ranks of each pop- 
ulation and we compare the sums. Again, if they are very 
different, the hypothesis that the two populations are drawn 
from the same parent distribution is ruled out. This test 
strongly supports the hypothesis that early-types and ellip- 
ticals have a different median at different redshifts, while for 
SOs it cannnot exclude the similarity of them (giving, respec- 
tively, a probability of 0.06%, 1.17% and 17.12%). We note 
that these results are fully in agreement with the bootstrap 
errors (see Table [Sjl. 

A detailed summary and comparison of the results of 
the different tests is shown in Table \KT\ both for the mass- 
limited sample and the magnitude- "delimited" ones (see Q. 



APPENDIX B: HOW CAN RESULTS FROM 
DIFFERENT SAMPLES BE RECONCILED? 

From §[3]and §3] we draw different results. Summarizing, we 
have found that in our mass- limited sample there is no clear 
trend between ellipticity and mass for SOs. For ellipticals, 
this trend is only hinted, with more massive galaxies having 
slightly lower values of ellipticity, while for early-types it is 
quite strong and mostly due to the fact that ellipticals and 
SOs are found in different proportions at different masses. 

Comparing the ellipticity distributions at the two red- 
shifts, we have found an evolution for the early-types, with 
WINGS galaxies being proportionally more flattened than 
EDisCS galaxies. No strong evolution has been detected for 
ellipticals and SOs separately, except for the likely presence 
of an enhanced population of round ellipticals at low-z. Note 
that this trend for round ellipticals is opposite to the trend 



To take into account the WINGS incompleteness, here we 
consider rounded WINGS' weights, so that WINGS galaxies can 
weigh 1, 2 or 3, according to their real weight. 
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test 




early-types 






ellipticals 






SOs 




KOLMOGOROV- 


MASS; 


0.04% 




MASS: 


5.84% 




MASS: 


>20% 




SMIRNOV 


MAG: 


7.7% 




MAG: 


0.41% 




MAG: 


11.9% 








< 1% in 60% 






< 1% in 3% 






< 1% in 2% 






MASS: 


1-5% in 25% 
5-10% in 7% 




MASS: 


1-5% in 11% 
5-10% in 12% 




MASS: 


1-5% in 12% 
5-10% in 10% 




MOSES 




> 10% in 8% 






> 10% in 74% 






> 10% in 76% 




simulations 




< 1% in 28% 






< 1% in 16% 






< 1% in 1% 






MAG: 


1-5% in 25% 
5-10% in 17% 
> 10% in 30% 




MAG: 


1-5% in 23% 
5-10% in 10% 
> 10% in 51% 




MAG: 


1-5% in 9% 
5-10% in 13% 
> 10% in 77% 




MANN-WHITNEY 


MASS: 
MAG: 


0.06% 
13.18% 




MASS: 
MAG: 


1.17% 
0.07% 




MASS: 
MAG: 


17.12% 
3.4% 


7^ 



Table Al. Summary of the results of different tests applied both to our mass and magnitude- "delimited" samples (see ^4}. The K-S, 
being as general as possible, gives an indication on whether two distributions can derive from the same parent distribution; the Moses 
tests the equality of scale parameters when the assumption of common medians is not reasonable (hence it tests the shape); the Mann- 
Whitney tests if there is a shift in median of the two population (hence it tests the compatibility of the medians) . For the mass-limited 
sample, only galaxies above log AZ/Mq 10.2 are considered. Moreover, for WINGS we also take into account the completeness weights 
(see text for details). The symbol 7^ means that the considered test can state that the two populations are drawn from different parent 
distribution, while the symbol = means that the considered test is not conclusive. 



for all early-type galaxies (rounder vs flatter at low-z, respec- 
tively), therefore we must be observing two distinct evolu- 
tionary effects. 

In contrast, from the analysis of the magnitude- 
"delimited" sample, we cannot exclude that, both in the 
case of early-type galaxies and of SOs, the galaxy samples 
at high- and low-z are drawn from the same parent distri- 
bution, although the change in the median ellipticity values 
with time seems to indicate an evolution instead. Moreover, 
we have found a significant evolution {2a error) of the el- 
lipticity distribution of elliptical galaxies, due mainly to a 
different median of the distributions. 

To understand the reasons for the observed discrepan- 
cies, we have compared directly the ellipticity distributions 
of the mass- limited and magnitude- "delimited" samples for 
the same type of galaxies at the same redshift (plots not 
shown) . 

The K-S suggests different distributions {Pk~s ~ 0%) 
for WINGS early-types and SOs, while it is inconclusive in 
all other cases (i.e. WINGS ellipticals; EDisCS early-types, 
ellipticals and SOs, Pk-s » 20%). Going deeper into the 
analysis, WINGS early-types show incompatible values both 
of median (the Mann- Whitney test gives a probability of 
0.90%) and of scale parameter (the Moses test gives a prob- 
ability < 5% in 80% of the simulations respectively), while 
WINGS SOs have different scale parameters (the Moses test 
gives a probability < 5% in 83% of the simulations respec- 
tively). 

The origin of the observed differences in the distribu- 
tions probably lies in the fact that galaxies in the two sam- 
ples are characterized by different properties; in particular in 
the magnitude- "delimited" samples, selecting galaxies only 
in the magnitude range —19.3 > Mb -1-1.2082 > —21, we are 
loosing (the few) most massive galaxies, and a large fraction 
of the less massive galaxies. 

This is evident in Figure IBll where we compare the 
galaxy stellar mass functions of our mass- and magnit ude 
(de)limited samples, as derived in I Vulcani et all (|2010bt l. It 



is clear that the magnitude- "delimited" sample is incomplete 
at low and at high masses, and much more so for SOs at low 
masses in the local Universe than at high-z. 

The net effect of the differential incompleteness in mass 
at high- and low-z in magnitude selected samples is to arti- 
ficially deprive the low-z distribution preferentially of high 
ellipticity (SO) galaxies. The loss of low mass SOs (more flat- 
tened) at low-z greatly reduces the differences between the 
high- and low-z ellipticity distribution of early-type galax- 
ies, bringing their medians to be almost consistent and the 
K-S test to be inconclusive. 

For ellipticals, the net effect of the incompleteness of 
the magnitude- "delimited" sample is to exacerbate the dif- 
ferences with redshift, again subtracting low-mass (hence 
higher ellipticity) ellipticals at low-z. 

The incompleteness in the mass distributions of the 
magnitude- "delimited" sample therefore seems to be con- 
sistent with the differences we observe between the elliptic- 
ity distributions of the mass-limited and of the magnitude- 
"delimited" sample. The magnitude- "delimited" sample is 
biassed; in particular, early-type galaxies on the red se- 
quence and with —19.3 > Mb + 1.208z > —21 are not 
representative of the overall population. 
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